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Abstract 

The first-order phase transition of 0(3) symmetric models is considered in the hmit of 
high temperature. It is shown that this model supports a new bubble solution where the 
global monopole is formed at the center of the bubble in addition to the ordinary 0(3) bub- 
ble. Though the free energy of it is larger than that of normal bubble, the production rate 
can considerably be large at high temperatures. 



Since a general theory of the decay of the metastable phase was developed [jl], H, the study 
of first-order phase transitions have attracted the attention due to their possible relevance 
to the physics of the early universe. The semiclassical expression of the tunneling rate from 
false vacuum to true vacuum is given by the bubble solution of lowest (Euclidean) action 
both for the first-order phase transitions at zero temperature field theories P| and those at 
finite temperature 0, However, the above case does not include the possibility that there 
exist various decay modes between two given classical vacua. In this note, we explore such 
problem by examining a scalar model of internal 0(3) symmetry at finite temperature and 
examine a new bubble solution with a global monopole. 

Suppose that the model of our interest contains a series of stationary points which support 
bubble solutions, a decay probability per unit time per unit volume is given by 

^/V = ^A„e-^"/^ (1) 

n 

where n represents the n-th local minimum. For each n, -B„ is given by a value of Euclidean 
action for n-th bubble solution and An is estimated by integrating out the fiuctuations around 
a given n-th configuration. 

We choose the 0(3)-symmetric scalar models described in terms of an isovector field 
(f)"' {a = 1,2,3) as the first possibility for sample calculations. At finite temperature the 
Euclidean action is 

s.^[dt.j<e. + vir)^, (2) 

where (3 = h/ksT. Although our argument followed is general and does scarcely depend on 
the detailed form of scalar potential if it has the vacuum structure with both a true and a 
false vacuum, we will consider a specific model, i.e. a 0^-potential 

V(0) = 4(0' + - v'f + Vo (3) 
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where denotes the amphtude of scalar fields 0" defined by = y/ (/)"■(/)"■. Here we consider 
the transition from the symmetric vacuum to the broken vacuum, i.e. < a < 1/2, and 
choose Vo as — Aaf^ in order to make the value of 5*^; coinside with Bn in Eq.(|l]). This does 
not lose generality since the case of the transition from the broken vacuum to the symmetric 
one (—1 < a < 0) also possesses the same type of bubble solutions if we replace to f — 0. 

It has been known that the rate of vacuum transitions at finite temperature has relevant 
amount of contributions from 0(3)-symmetric sphaleron-type bubbles [§], In high tem- 
perature limit this contribution dominates and then we can neglect the dependence along 
tg-axis. Here we are interested in the regime where thermal fiuctuations are much larger 
than quantum fiuctuations and address the problem for different kind of finite temperature 
bubbles which connect global 0(3) internal symmetry to that of spatial rotation. We now 
ask for a solution of the field equations that is time-independent and spherically symmetric, 
apart from the angle dependence due to the mapping between the (6*, ip) angles in space and 
those of isovector space 05^(= 0"/0) such as 

05^ = (sin n6' cos ny), sin n9 sin rup, cos n9), (4) 

where the allowed n is or 1 in order to render the scalar amplitude a function of r only. 
Under this assumption, the Euler-Lagrange equation becomes 

fl + l^-S r-(j) = — (5) 
dr'^ r dr c?0 ' 

where 6ni in the third term of Eq.(|^) denotes Kronecker delta for ra = or 1. 

The condition that the theory should be in false vacuum at spatial infinity fixes the 
boundary value of field, i.e. limr^oo 0. To be nonsingular solution at the origin of 
coordinates, the boundary condition is 

[ 0(0) = if n = 1. 



When n = 0, it is well-known bubble solution at finite temperature We will analyze 
n = 1 case and show that there always he n = 1 solution if the equation contains n = 
solution. A brief argument of the existence of n = 1 solution is as follows. 

If we regard the radius r as time and the scalar amplitude 0(r) as the coordinate of 
a particle, Eq.(|D describes a one-dimensional motion of a unit-mass hypothetical particle 
under the conserved force due to the potential —V{(f)) and two nonconservative forces, i.e. one 
is the friction of time-dependent coefficient — f ^ and the other is time-dependent repulsion 
^0. Hence, in the terminology of Newton equation, n = 1 solution in Fig. 1 is interpreted 
as the motion of a particle that starts at time zero at the origin (0(0) = 0), turns at an 
appropriate nonzero position at a certain time {4>(tturn) = (pturn), and stops at the origin at 
infinite time {4>{oo) = 0). 

At first, if one considers a hypothetical particle at the origin at time zero, it is accelerated 
by the time-dependent repulsion of which coefficient is divergently large for small r. Since 
one initial condition is fixed by the starting point (0(0) = 0), such motions near the origin 
are characterized by another initial condition C 



where C is the initial velocity of a particle which should be tuned by the proper boundary 
condition at infinite time (0(oo) = 0). From now on let us consider a set of solutions specified 
by a real parameter C and show that there always exists the unique motion of 0(oo) = for 
an appropriate C. When C is sufficiently large, the acceleration near the origin due to time- 
dependent repulsion and conservative force is too strong and then the particle overshoots the 
top of potential —V{(f) = v) and goes to infinity (0(oo) = oo), despite the deceleration due 
to the friction. Since the solution of our interest is the motion which includes a return, we 
will look at the motions of which C is smaller than a critical value Ctop- Here Ctop gives the 
motion that the particle stops at the hilltop of the potential —V at infinite time 0(c>o) = v. 
Suppose that C be too small, i.e. C is smaller than another critical value Co (Co < Ctop), 




(7) 



3 



particle turns at a point smaller than 0o (0 < 0o < v) where — V^(0 = 0o) = and thereby 
it can not return to the origin. Moreover, the particle which turns at a point too close 
to V arrives at the origin at a finite time, since the time- dependent friction and repulsion 
during returning do not play a role due to the vanishing of their coefficients 1/r and 
while the particle stays near v for sufficiently long time. By continuity for an appropriate 
C (Co < C < Ctop) there is a solution which describes the motion that starts at time zero 
at the origin turns at a position (ptum between 00 and f , and then stops at the origin at 
infinite time. Here for the existence proof oi n = 1 solutions, we have used the properties 
of scalar potential that it has a local minimum at = 0, a local maximum at (pbottom 
(0 < (pbottom < 0o) and an absolute minimum at v. It is exactly the same condition as that for 
the existence of n = solution. Hence we completes our argument that the equation always 
contains n = 1 solutions when it has n = solution. A rigorous proof for the existence of 
n = 1 solution is demonstrated in ref . . We have no analytic form of solutions, so numerical 
solutions for n = and n = 1 bubbles are given in Fig. 1. 

Suppose that each classical solution describes a critical bubble corresponding to a different 
decay channel, from the profile of energy density Tq (see Fig. 2) we can read the following 
characteristics of n = 1 bubbles. First, the boundary value of scalar field at the origin 
always has that of false vacuum, i.e. 0(0) = 0, so this implies that there remains a false- 
vacuum core inside the true- vacuum region of the bubble due to the winding between 0(3) 
internal symmetry and spatial rotation. Second, the value of energy density has a local 
minimum at the origin, i.e. 3C^/2 + 1^(0) in which C is a constant appeared in Eq.(|^), 
increases to the maximum at Rm in Fig. 2 and decays below V{0). This implies that a 
matter droplet is created inside the bubble and is surrounded by inner bubble wall with size 



of order Rm ~ ^/rriHiggs = \/4(3 + 2a) Af. Third, if we read the expression of energy density 
when the scalar amplitude has maximum value = (ptum, its derivative term vanishes and 
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then it becomes 

T'o = ^ + V{<l>turn). (8) 

The potential term V{(f)turn) can be neglected in thin-wall limit, so the object inside the 
bubble has a long-range hair which penetrates the inner bubble wall. Since the central 
region of this matter is in false symmetric vacuum due to the hedgehog ansatz of scalar field 
in Eq.(D and for large r the scalar field goes to the true broken vacuum but the long-range 
tail of energy proportional to r has a cutoff at the outer bubble wall at Rn=i in Fig. 2, we can 
interpret the matter aggregate inside the n = 1 bubble as a global monopole of size Rm [^]. 
Even though there is a no-go theorem for static scalar objects in spacetime dimensions more 
than two ^ , the global monopole of this case can be supported inside the Euclidean bubble 
configuration as a smooth finite-energy configuration due to a natural cutoff introduced by 
the outer bubble wall. Fig. 2 shows that the radius of = 1 bubble is larger than that of 
n = bubble. It can be easily understood by the conservation of energy, i.e. the additional 
energy used to make a matter aggregate is equal to the loss of energy due to the increase of 
the radius of bubble. 

Though we do not have any analytic solution, we can estimate i?„ in Eq. (|^) by use of the 
obtained bubble configurations through numerical analysis for given parameters of theory. It 
has already been proved that n = solution describes the nontrivial solution of lowest action 
|§] and from Fig. 3 we read that the value of action for n = 1 solution, vB[/T, is larger than 
that of = 0, vB'q/T, irrespective of the shape of scalar potential. Fig. 3 also shows that 
the ratio B[/ B'q becomes small while the difference B[ — Bq increases as the size of bubble 
becomes large in comparison with the mass scale of theory, i.e. thin-wall limit. This can 
be understood as follows; the amount of energy consumed to support a global monopole at 
the center of bubble is proportional to the radius of bubble due to its long-range tail (see 
Eq-®), however the free energy to support the bubble itself, vB'q or vB[, is proportional to 
the cubic of bubble radius in thin- wall limit. 
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The next task is to estimate the pre-exponential factor An in Eq.(||). The scheme of 
computing Aq for n = bubble solution which is the lowest action solution with one negative 
mode was given in the second paper of Ref.0. Here, let us attempt to calculate the pre- 
exponential factor Ai for n = 1 bubble. Considering the small fluctuation around n = 1 
bubble S(j)n=i = J2k CfcV'fc, we obtain a Schrodinger-type equation for three particles in three 
dimensions 



IdV 



^\ = Afc^^ (9) 



It looks too difficult to solve this equation since it includes 6 and dependent off-diagonal 
terms in its Hamiltonian and the potential form is only determined numerically, however the 
eigenfunctions of six zero modes due to three translations and three rotations are explicitly 
given in a form 

and 

rk,i = Nre,,kx'V'rn=, = iV,e,,,fV„=i, (11) 

where Nt and Nr are normalization constants. Here we give few comments on the number of 
negative modes. First, every component of translational zero-mode eigenfunction in Eq . (|T0|) 



has {9, y9)-dependence where its radial part, — has single node at the origin, and 
each a = i component includes an additional {6, v9)-independent part, which has no 

node. Every a i component of rotational ones, 0n=i, has one node at the origin. For the 
sake of simplicity, let us examine the problem for the perturbation with specific direction; (i) 
one for amplitude (50^=i = r"'S(f)n=i) and (ii) two for transverse (50^=1 = {6"''' — 'r°"r^)6(j)n=i) ■ 
Under these fluctuations Eq.(P) reduces to 

(_V2 + U{rm = ^krk, (12) 



where U{r) = ^ for the first case (i) and U{r) = 

When f/(r) = , each single 'a' component of Eq.(|l^ is nothing but the fiuctuation 
equation for n = bubble which contains a nodless s-wave mode as the unique negative mode 
0, 1^. For the fiuctuation of scalar amplitude for n = 1 bubble, the lowest mode can not be 
nodless s-wave mode but / = 1 mode with single node at r = since ip'^ should have 6 and ip 
dependence proportional to f°- even though the operator in Eq. ([T^) is a scalar operator. It is 
analogous for the perturbation to the transversal directions. It implies that the eigenfunction 
of negative mode may take a somewhat complicated form which depends on angels. Second, 
the n = 1 solution which is a parity-odd bounce solution, 0(r = 0) = 0(r = oo), is supported 
by assuming the winding between three spatial rotations and those in internal space. It seems 
that the argument for the system of quantum mechanics with one time variable in Ref.0 
does not forbid directly the existence of n = 1 solution as a bubble configuration. We know 
very little on the counting of negative modes, and then there should be further work on 
this issue. Since the operator in Eq.(P) is even under parity transformation and is covariant 
under the rotation, i.e. — > 0'''^x^ and ^J^i^') — 0"'^il)\{x) where both O*-' and O"* are the 
elements of 0(3) group, vector spherical harmonics is a method to investigate the modes of 



for the second case iii). 



Pn = l 



which the eigenfunctions can be chosen to be either even or odd in parity |T0 . 

In order to compute the decay rate accurately, we should calculate the nonzero modes 
Afc, which is extremely difficult even for n = bubbles. However, the dimensional estimate 
may reach a rough result such as 



At high temperature limit, T >> v., both decay rates for n = and 1 bubbles, of course, 
increase and, moreover, the relative decay rate of n = 1 bubble to that of n = is also 
enhanced exponentially. When v/T ~ 10^^, the order of T^^"* /V^^^ is around 0.17 in a 
thin-wall case (A = 1, a = 0.12) and it is around 28 in a thick-wall case (A = 1, a = 
0.47). Therefore, at high temperatures and in thick-wall case, this n = 1 bubbles can be 



preferred to those of n = 0, and obviously the existence of another decay channel can enhance 
considerably the total nucleation rate of bubbles, T = T^^^ + T^^\ except for the bubbles with 
extremely thin wall. 

The main consideration of the paper was to gain an understanding as to how bubbles 
with global monopoles nucleated in first-order phase transitions at high temperatures and 
the next question will be how the high-temperature bubbles grow particularly at the 
site of global monopole. Once we assume these bubbles in early universe where the gravity 



effect should be included, the bubbles with solitons may result in interesting phenomena [|T2| 
in relation with inflationary models, i.e. the inflation in the cores of topological defects 
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Figure Captions 



FIG. 1. Plot of bubble solutions, n — and n — 1 configurations are shown as dotted and 
solid lines, respectively. The parameters chosen in the figures are: A = 1, a = 0.12, and 
Vo = -0.12^;^ 

FIG. 2. Plot of energy density Tq . The parameters chosen in the figures are: A = 1, a = 0.12, 
and Vo = —0.12^"^ which is the minimum of energy density. 

FIG. 3. Plot of action Se (or equivalently Bn) as a function of a. Another parameters are 
chosen as A = 1 and Vq = —Xav"^. 
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Figure 1 
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Figure 2 
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